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Ginzburg-Landau theory of the zig-zag transition in quasi-one-dimensional classical 

Wigner crystals 
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Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020, Antwerpen, Belgium 

We present a mean-field description of the zig-zag phase transition of a quasi-one-dimensional 
system of strongly interacting particles, with interaction potential r~ n e~ r ^ x , that are confined by a 
power-law potential (y a ). The parameters of the resulting one-dimensional Ginzburg-Landau theory 
are determined analytically for different values of a and n. Close to the transition point for the 
zig-zag phase transition, the scaling behavior of the order parameter is determined. For a = 2 the 
zig-zag transition from a single to a double chain is of second order, while for a > 2 the one chain 
configuration is always unstable and for a < 2 the one chain ordered state becomes unstable at a 
certain critical density resulting in jumps of single particles out of the chain. 
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I. INTRODUCTION 



During the last two decades the interest in self- 
organized systems has increased enormously both ex- 
perimentally and theoretically, due to its importance in 
solid-state physics, plasma physics as well as in atomic 
physics. Wigner crystals are an elementary example of 
self-organization which has been realized in very diverse 
systems as e.g. electrons on liquid heliurol, by using the 
well-known Paul and Penning traps^ to confine ions in 
a limited region, in dusty plasma^, and more recently us- 
ing static and radio- frequency electromagnetic potentials 
where crystallization was realized through laser cooling^. 
Additionally it has been proposed that these structures 
can be used for a possible implementation of a scal- 
able quantum information processor^, and as quasi-one- 
dimensional (Q1D) Wigner crystals^. The theoretical 
analysis of these crystal structures has been realized pre- 



viously, for 3D&±ii, 2EUirii and Ql 



systems. From 



those studies it was shown that structural phase transi- 
tions can be induced by varying the strength of the exter- 
nal confinement potential and/or the density of particles. 

In the present paper we concentrate on the ordered 
state of identical particles that are confined in a Q1D 
channel. When the particles move in 2D and are con- 
fined by a parabolic^ or hard walU& potential in one 
of the in-plane directions the particles arrange them- 
selves in parallel chains at low temperature. Previously, 
it was founcU^ that with increasing density (or decreas- 
ing strength of the confinement potential) the system 
passes through a sequence of first and second order phase 
transitions where at each point the number of chains 
changes. Of particular interest to us is the one chain to 
two chain transition which for a parabolic confinement 
potential was found to be a second order phase transi- 
tion which occurs as a zig-zag transition. Such a tran- 
sition was observed cxperimentallyi£~— in systems with 
a finite number of particles, and the effects of a narrow 
channel and finite size of the system on the diffusion was 
recently analyzed in Refs. [2l[ and [22j|. Recently it was 
found theoretically^ that the analytic form of the con- 
finement potential is very important for the occurrence 



of the zig-zag transition and the order of the phase tran- 
sition. Therefore, in the present paper we generalize the 
previous analysis to an arbitrary power law confinement 
potential (i.e. y a ) and also to arbitrary inter-particle in- 
teraction which we model by r~"e _r / A , which simulates 
most of the relevant experimental particle-particle inter- 
actions. With this model potential we can simulate both 
short range and long range interactions. 

We are interested in the behavior of the system at the 
zig-zag transition i.e. at the critical point. This can 
be viewed as a spontaneous symmetry breaking and we 
will cast the problem into a mean field theory based on 
Landau's theory of phase transitions. In this way we 
will construct a Ginzburg-Landau theory, for the single 
to two chains transition in a quasi-one-dimensional sys- 
tem of interacting particles. We generalize the approach 
of Refs. [lH and [24| to arbitrary power-law confine- 
ment and inter-particle interaction potential. We ob- 
tain a Ginzburg-Landau equation for the order param- 
eter close to the transition point, and determinate all the 
relevant parameters in this equation. The order param- 
eter is the distance of the particles from the trap axis. 
By considering a large number of particles and using the 
local density approximation, we can consider the crystal 
as a continuum, so that the order parameter becomes a 
field. 

The present paper is organized as follows. In Sect. [TT] 
we describe the model system and using Landau theory 
we find the behavior of the system close to the transition 
point. Next we derive a Ginzburg-Landau equation for 
the system finding the dispersion relation. In Sect. IIIII 
the results for the critical point and the normal mode 
spectrum are discussed. Our conclusions and a discussion 
of possible quantum effects are given in Sect. HVl 



II. THEORETICAL FRAMEWORK 

We consider a two-dimensional system consisting of N 
particles with mass m and charge q, which arc allowed to 
move in the x — y plane. The charged particles interact 
through a repulsive interaction potential; they are free 
to move in the x direction but are confined by a one- 



2 



dimensional potential which limits their motion in the y 
direction. The total energy of the system is given by: 



1 N 

E = y m 22 *i + Vcon f + Vint > 



(1) 



where V con f and Vi n * are the confinement and interaction 
potential respectively, given by 



R c 

N N 
y ' ^ ] ^pair ( r ij , 

2—1 J>2 



(2a) 



(2b) 



where V pa i r (r^ ) represents the inter-particle interaction. 
The latter one will be taken as a screened power-law po- 
tential as follows: 



r, 



q 2 R r 



pair 



eR 



(3) 



where r»j = |r, 



Tj l represents the relative position be- 



tween the z-th and the j-th particle, the exponent n is an 
integer and e is the dielectric constant of the medium the 
particles are moving in. In the above, R is an arbitrary 
length parameter which we introduced to guarantee the 
right units. The energy can be written in dimensionless 
form 



N 
i=l 



N N 

EE 

i— 1 j>i 



(4) 



with dimensionless frequency v given by v = 
vt/ujQ, while ujq measures the strength of the con- 
finement potential and to = 1/wo is the unit 
of time. The energy is expressed in units of 

E Q = (m Wo 2 /2) Tl/(n+Q) ( g 2 /e) Q/(Tl+Q) J R( 2 "- a )/(»+«) and 
all distances are expressed in units of tq — 
(2< Z 2 /mw 2 e) 1/( " +Q) i?("+ Q - 3 )/("+ Q ). Additionally, the 
dimensionless parameter k = tq/X represents the screen- 
ing parameter of the potential. Limiting cases of this in- 
teraction potential are: Yukawa potential (n = 1), power- 
law potential {k = 0), Coulomb potential (k = 0, n = 1) 
and dipole interaction («; = 0, n = 3). We introduce a 
dimensionless linear density r\ defined as the number of 
particles per unit of length along the unconfined direc- 
tion. 

In Ref. [23| it was demonstrated that for a > 2 the 
one-chain configuration is not stable for any values of 77 
and v. Only for a = 2 the system exhibits a continuous 
transition from the one-chain to the two-chain configura- 
tion (i.e. zig-zag transition) at a transition point defined 
by a critical density (rj c ) or a critical frequency (v c ). For 
a < 2 the ground state configuration of the particles is, 
below rj c or above v c , arranged in a single chain. Beyond 
this critical point particles are expelled one by one from 
this chain to positions parallel to the chain. 



For the special case of a = 2 and before the transi- 
tion point, the particles crystallize around the minimum 
point of the confinement potential V , at the positions 
r imear _ (j/ 7 y)e x with i an integer. The stability of the 
linear chain along the x axis requires a relative trans- 
verse trap frequency exceeding a threshold value v c or a 
linear density smaller than rj c . At this critical point, the 
configuration has a structural instability, such that for 
v < v c or 77 > j] c the particles are organized in a zig-zag 
structure, ordered in two chains with equilibrium posi- 
tions Y z i igza9 = (i/r])e x + (~l) l (c/r])e y , where c is a real 
and positive constant, d = c/rj represents the distance of 
each particle from the confinement potential minimum 
and D = 2d indicates the lateral separation between the 
two chains. 



A. Landau Theory for the zig-zag transition 

For the case a = 2 and near the zig-zag regime we 
follow the Landau theory approach of Ref. (23(, and 
expand the total potential energy of the system as a 
function of the order parameter c in a polynomial i.e. 
V(c) = Vi c h — Ac 2 + Be 4 , where V\ c h represents the po- 
tential energy for the one-chain configuration. Minimiz- 
ing the potential energy we obtain the condition 

- K 2l±± 2 

2£ [mf + K V n ~ 1 (2j + 1)] — " - V - = 0. (5) 

From this equation we obtain the value of r) c (v) or 
v c (rf) at which the single chain configuration becomes 
unstable. Considering r\ c and expanding the potential 
energy around this critical value, we find that the order 
parameter close to the transition point is given by 



c = Y(v : n, K)\r) — r] c \'' 



(6) 



with Y(v; n, k) 



VY B Y C - Y A Y D /Y B where Y A 
2 /o^2 v „ _ 2n ( n + 2)r]^S 4 



V^kSi - u 2 /2?7 2 , Y B 



2(2n+l)r 1 ^ 1 K S 3 + 2r 1 ^- 2 K 2 S 2 , Y c = n 2 i 1 ^- 1 S 2 + (2n - 
I^-^Si+^-Wo+t; 2 /^, Y d = 2n 2 {n+2)i 1 ™- l S A + 
2[n(3n+l)-l]?7™- 2 KS , 3+2(3n-l)77™- 3 K 2 S , 2-l-2?7™- 4 K 3 1 S*i, 
with 



E 



(2i + l)"+ fe ' 



(7) 



The value of Y = Y(v; n, k) is plotted in Fig. [T] as 
function of k for different values of n and a fixed value 
v = 1. Notice that the curves Y(n, k) for different values 
of n cross each other at some value of k. The critical 
exponent of the order parameter of the zig-zag transition, 
Eq. ([B]), was verified experimentally on a low dimensional 
dusty plasma in Ref. 
fil, El and M. 



251 ] and theoretically in Refs. [2c 
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FIG. 1. (Color online) The proportionality coefficient 
Y(v; n, n) as a function of n for different values of n with 
v = l 



B. Ginzburg-Landau Lagrangian for the zig-zag 
phase transition 

Recently, this zig-zag transition (for a = 2, k = 
and n = 1) was cast into a mean- field description re- 
sulting in similar expressions as in the Landau theory 
of phase transitions. This resulted in a one-dimensional 
Ginzburg-Landau type non-linear field theory 2 ^. Here, 
we will extend the previous calculation to the more gen- 
eral problem described by the energy Eq. (j4|). We start 
by considering the system in the situation that the one- 
chain configuration is stable but that it is close to the 
transition point. The equilibrium positions of all the par- 
ticles are along the x-axis. We consider small oscillations 
around the equilibrium position of each particle as follows 



(i/v) 



and y" 1 = + i/j. Then the relative 



position between the particles can be written as follows 

1 1/2 



1 



A, 



(^■) 1/2 , 



(8) 



with Aij = (Ay) 2 , Tij = 2Aij(xi - Xj) and ey = (x t - 
Xj) 2 + {yi — yj) 2 , where Ay = (i — j)/r). Now, we assume 
that the vibration amplitudes in the axial and transverse 
direction are much smaller than the distance between the 
particles, i.e. ry, ey <C Ay. We expand Eq. and the 
exponential term as follows 



r y ~ Ay ~\~ {xi Xj ) 

(xi - Xj) (Xi - Xj) 2 \ (yt - y-j) 2 



1 - 



(Ay) 2 



2A, ; 



(Xi Xj) 



(9) 



. {Vi-Vj? 
2A„- 



(10) 



and similar the n-th power of the inverse of Eq. 
i.e. 1/V™ , as a Newton binomial around the equilibrium 
positions. These expansions result in a decomposition of 
the total potential as 



int i 



(ii) 



where the label I indicates the order of the expansion. 
Each order of the expansion of the interaction potential 
can be written as 



V, 



OO 



(0 
y ' 



(12) 



where the expansion terms up to fourth order are given 
by 



4 



(0) 



'J 
(1) 



(2) 



5 = 



(4) 



-«Aj 



(A y - 

e -«A„ 

2(A~ 



*J 

y^[$ijn + KAij}(Xi -Xj), 



[{n{n + 1) - 2^ ij nKA ij + k 2 A|) (x* - Xj) 2 - (n + kA^ ) (j/j - 2/j) 2 ] , 
y^n(x t - x,) |- {-f(n 2 + 3n + 2) + (n + 1)kA 13 + W,,/rA 2 ,J (a* - x,-)" 



2(A 2J 

+ (i%(n + 2) + fcA,j ) (y< - yj -) 2 
1 



2(A ij )"+ 4 



12 



{n 6 + 6n^ + lln + 6) + 3(n 2 + 3n + 2)kA ?;j + 2(n + 1)k A£ (a* - Xj-) 



-(n 2 + 5?i + 6) + (9(n + 2)0y + 2{n + l))/sAy + 2k 2 A 2 - ) (x* - ay) a (i/j - Vjf 



+ Q(n + 2) + - %) 4 



r 



With -dij = Ay / 1 Ay | . It is sufficient to restrict our- 
selves to terms up to the fourth order and thus the po- 
tential can be written as V int « Vf£} + v£? + V^) + 



V; 



(3) 



V; 



(4) 



Representation in reciprocal space 



Now we assume that the particles are pinned in 
the longitudinal direction and that they can only os- 
cillate in the transverse direction (xj = 0). There- 
fore, their normal axial modes can be neglected 
and wc can discard the coupling to the longitudi- 
nal modes. In this regime, we find that W/, 0) = 



e-* A -7(A"), W jf = 0, Wi 

y j re-^/(2(A ij r+ 2 ),W i 
n [(n + 2) + 8/tAy] ( W - y 3 ) 



(2) 



(n+ kAjj) (yi 



(3) 



= 0, and W, 



(4) 



4_-kA. 



V(8(A lJ )" +4 ) 



In order to find the representation in reciprocal space, 
we define the normal modes of vibration in the transver- 



sal direction with wavevector fc as = ip 



(+) 



with amplitude 



;(+)2 



Following the 



standard process to find this representation as shown in 
e.g. Rcf. and using Plancherel's theorem 2 ^ for the 
confinement potential transformation, the different terms 



of the potential become 

,,2 



conf 



(1) 

n t 



V; 



(3) 



V; 



int 
(4) 



fe>0 

0, 

k>0 

0. 



(13a) 

(13b) 
(13c) 

(13d) 



■ill I 



£ ^1^2,^3,^4) IJ ^ m (13e) 

fei+fe2+fc3+A;4=0 m=l 

where aij_(n, «, fc) 2 = — w(n, k, fc) 2 , A(fci, fe, fca, foi) = 
n77 n+4 [(n + 2)A (k 1 ,k 2 , fc 3 , fc 4 ) + 8RA 1 {k 1 ,k 2 , k 3 , **)]/8, 
with ii ~ k/tj and 



k, fc) 2 = 2 77 



n+2 



e-J" . 2 / fc 
Z^-^T2- sm I Jo 



CO 

3=1 



— Sill 7 — 
-1 ^ J 9 



,(14a) 



A (h,k 2 ,k 3 ,k 4 ) = n sin(i^), 



iV ' — ' 7' 

j>0 J 



m— 1 
4 



A 1 (fc 1 ,fc 2 ,fc 3 ,fc 4 )=^EliT3 ri sin u^)' 



j>0 J 



(14b) 



Due to the condition that the motion of the particles 
are restricted to the longitudinal direction it becomes ap- 
parent that the first and third order term of the interac- 
tion potential will be zero, and additionally we know that 
the first derivative equals zero because it is the necessary 
condition to have an equilibrium configuration. 



5 



2. Minimum frequency of the interaction potential 



around this value (fc = fco — Sk), and wc obtain 

ui±(n, «, fco — 8k) 2 = uj±(n 7 k, fco) 2 

+h(n,k) 2 Sk 2 + 0{Sk 4 ), (15a) 

A(k u fc 2 , fc 3 , fc 4 ) = ^Mn, ft) + 0(Sk 2 ), (15b) 

From the definition of uj±(n, ft, k) 2 we find that its min- 
imum value is located at fc = tt. Lets expand for fc where uj±(n, k, fc ) 2 =— zu(n, k) 2 with 



,^2 



U7(n, k) 



n+2 



2n$ I e~ 2f \ n + 2 , - I + 4k$ ( e~ ZK , n + 1 



h(n,R?={%) ^(-l)^ +1 e-^ 

3=1 



-2k „ 



.4(n, k) 



n+4 



-$ e _ ^ K ,n, - + k$[ e~ iK ,n- 1, - 



1 



- + 1 ) $ ( e~ ZK , n + 4, - I + 8k$ [ e~ 2 ~ K , n + 3,l 



(16a) 
(16b) 
(16c) 



where <&(z,s,a) is the Lerch transcendent defined as 
$(z,s,a) = J2kLo zk I ( k + a Y '• In Table U we show the 
limiting behavior of these terms. It is important to note 
that the square of the transverse frequency is negative. 



3. Stability of the system 

The system is stable when the second order term of 
the total potential energy (i.e. the coefficient of if) 2 .) is 
minimum. For the parabolic case (a = 2) the confine- 
ment potential term contributes to the second order of 
the total potential energy, then 



y( 2 >=]>> ± (n,ft,fc)V 



k>0 



where the transverse frequency is 

u>±(n, ft, fc) 2 — v 2 — w{n, k, fc) 2 



(17) 



(18) 



and we note that its minimum value is reached for fco = it. 
Thus the critical value for the confinement frequency is 
given by 



where 5 v (n, k) = v — v c (n,k) . In the limiting case 
of a Coulomb inter-particle potential (k = 0, n = 1) 
and considering 77 = 1, we find v c (l,0) = ^/7C(3)/2 = 
1.45038, h(l,0) = y/log{2)/2 = 0.58871 and -4(1,0) = 
93£(5)/64 = 1.50679, which agrees with the results of 
Refs. 0] and 0- In Table [H we show the values 
of these terms for different interaction potentials. No- 
tice that the critical confinement frequency v c decreasing 
with increasing screening ft, and it increases with increas- 
ing density. The relation between v c and n depends on 
the density value, as will be discussed later. 

Additionally, from a simple expansion of the disper- 
sion relation Eq. (|18[) around the equilibrium positions 
and for values of the frequency and density close to their 
critical values, the value of the parameter c can be found 
from the non-linear algebraic equation 



v 2 - 2n n+2 



L fe[(2j-l) 2 +c 2 ]- 



K h Pi-i) 2 + c 2 ]^ 



(21) 



v 2 (n, k) = w(n, k) 



(19) 



When v > v c the ground state configuration is a one- 
chain organization of the particles. For v < v c the linear 
chain is unstable and the particles are arranged in a two- 
chain structure through a zig-zag organization. When 
v is sufficiently close to the critical value v c , an effective 
potential can be derived for the transverse normal modes 
■0fe with wavevector fc = fco — Sk, such that 8k <C 1. The 
second order term of the effective potential is given by 
Eq. (fT7|) . where its coefficient, Eq. (|15a|) . can now be 
written as follows 

uj±(n, ft, fco — 8k) 2 = S v (n, k) + h(n, k) 2 5k 2 , (20) 



4- Continuum approximation 

Close to the transition point (dk -C 1), the trans- 
verse deviation of the particles is very small and we can 
use a continuum approach for these modes. In doing 
so we replace the discrete sum over fc by an integral 
Sfc — * / d(5k)N/2Tr. Using the Fourier transform we 
obtain a continuous form for the modes ip{ x ) as follows 
ipk = / dxip(x)e~ lSkx I VN . Then the remaining terms of 
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TABLE I. Behavior of the coefficients in the Ginzburg-Landau equation in two limiting cases, where k — k/t). 



vj(ii, ay 



h(n, k) 



A(n, k) 



ii < 1 



2r> 



E 



3=0 (3 + 1/2)"+^ 



f (i) n E^o(l-2iA) 



" I 2 3 1 M2^ ^j=0 (j+l/2)"+4 



K > 1 



2nr) rl 



E 



1 + 4 . 



i=0 (j + l/2)"-l 



E 



.3=0 (j + l/2)"+4 



TABLE II. Values of the critical parameters for different values of n and 









n = 0.50 






77 = 1.00 






77 = 1.50 




n 


K 


v c (n, k) 


h(n, k) 


A(n, k) 


v c (n, k) 


h(n, k) 


A(n, k) 


v c (n, k) 


h(n, k) 


A(n, k) 


1 


0.5 


0.43113 


0.38151 


0.06332 


1.36621 


0.57569 


2.13008 


4.02285 


0.82789 


62.6067 


I 


1.0 


0.31885 


0.30221 


0.04019 


1.21941 


0.53954 


2.02615 


3.86422 


0.81415 


68.1625 


I 


2.0 


0.15131 


0.15008 


0.01002 


0.90183 


0.42740 


1.28595 


3.44901 


0.76302 


64.8368 


2 


0.5 


0.37195 


0.35188 


0.06901 


1.74736 


0.80730 


4.85728 


7.52638 


1.71598 


299.492 


2 


1.0 


0.26019 


0.25394 


0.04230 


1.48781 


0.70376 


4.41649 


6.98946 


1.61460 


310.865 


2 


2.0 


0.11720 


0.11676 


0.01030 


1.04076 


0.50788 


2.70688 


5.95124 


1.40752 


282.655 


3 


0.5 


0.30339 


0.29542 


0.05605 


2.06259 


0.99201 


8.19081 


12.7470 


3.04303 


1047.17 


3 


1.0 


0.20566 


0.20332 


0.03331 


1.71625 


0.83557 


7.17465 


11.6677 


2.80582 


1048.42 


3 


2.0 


0.08952 


0.08936 


0.00794 


1.16342 


0.57508 


4.26315 


9.70860 


2.36334 


918.354 



the potential becomes 



V ^) = l f v 2 \^(x)\ a dx, 



V 



(2) 



—win, k) 2 tp{x) 2 



+h(n,n) [ — 



2-1 



dx, 



y [A) = o / A{n,K)ip{x) A dx. 



(22a) 

(22b) 
(22c) 



Finally, we obtain the Lagrangian L = J C(x)dx where 
the Lagrangian density C{x) reads 



= 5 [( d Mx)f - Kn, k) 2 {d x ^{x)f + w(n, k) 2 ^{x) 2 - v 2 \^{x)\ a - A{n, k)i>(x) A ] . 



(23) 



In the special case a = 2, k = 0, n = 1 this Lagrangian 
density is the one found in Ref. [24|], and it has the form 
of a Ginzburg-Landau equation (Refs. [24[ and [291). 
Defining tp(x) = rjip{x) and v(n, k) 2 = v(n, k) 2 jr\ 



n+Q 



we 



may find from Eq. (|23[) an expression for the potential 
energy density V: 



2V 

77™ 



K(n,k)ip(x) 4 + i(n,k) 2 \ip(x)\ a 
^(n, k)ip(x) 2 , 



(24) 



I 

and decrease with increasing k. Now, the density rj 
plays the role of a scaling parameter in the potential 
energy density (V), in the screening parameter (k), in 
the strength of the confinement (v) and in the order 
parameter (ip). For a = 2 we find the usual Landau 
energy expression for a second-order phase transition 

2V 

— = K(n, k)^{xf + [v(n, kf - Q(n, k)Mx) 2 . (25) 



with the real positive coefficients 0(n, k) = 
m(n,k) 2 /n n+2 and K{n,k) = A{n 7 k)/if+ 4 that 
are plotted in Fig. [2] as function of k for different 
values of n. Notice that both coefficients are positive 



From Eq. 
tp{x) as follows 



5. Equation of motion 
we obtain the equation of motion for 
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FIG. 2. Ginzburg-Landau coefficients for the potential energy density (Eq. Q24J)) as a function of the screening parameter 
k — k/tj, for different values of n. 



d?ip(x) - h(n, k) 2 d 2 x ip(x) 



cy.u 

w(n,k) tp(x) + —^-sign(ip(x))\ip(x)\ l 



+ 2A(n,k)i)(xf = 0. 



(26) 



In this context the order parameter ip{x) represents a 
continuous version for the value of c which is the dis- 
tance of the particles (d = c/rj) from the minimum of the 
confinement potential. When the order parameter varies 
slowly in space, the time independent version of Eq. (|26[) 
becomes 

2 

w{n, kfip - ^-sign^)]^- 1 - 2A{n, k)^ = 0. (27) 

We note that for a < 1 a one-chain configuration is 
not allowed, because ip = is not a solution of Eq. (|2T|) 
in this case. 

Considering a > 1 and defining a{n, k) = 
av 2 /2rj n+2 n(n, k) and b(n, k) = 2rj 2 K(n, k)/fl(n, k), the 
latter equation is reduced to 



1 — a(n, k)|-0| 



a-2 



b(n, k)ip 2 



0. 



(28) 



For a = 2 this equation results in a second order tran- 
sition, from the single chain (i.e. ip = 0) to the zig-zag 
(i.e. ip 0) configuration, with the critical point defined 
by a(n, k) = 1 which in fact is a generalization of Eq. 

For 1 < a < 2 and minimizing Eq. ([2"5]) we find 




FIG. 3. (Color online) Critical value of the exponent of the 
confinement frequency as a function of the density for different 
values of the strength of the confinement potential. We took 
the parameters n = 1 and k = 1 



2Q(n,k) 



2yyn+a c 



-K(n,k)j = 
(29) 

which represents a non-linear equation for the critical 
exponent of the confinement potential (a c ), which is the 
minimum value of a for which a one-chain configuration 
is the ground state configuration. From Eq. (f2T)|) we 
note that this critical value will be at most equal to 2, as 
shown in Fig. [3] for different values of the strength of the 
confinement frequency. 



Finally we may find analytical expressions for the order 
parameter from Eq. (|2"5)l for different values of a > 2, 
which are given in Table IIIII Notice that it is always 
possible to find a ip ^ which indicates that the single 
chain configuration is always unstable when a > 2. For 
a = 2 we find tp ^ only when a(n, k) < 1. 
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TABLE III. Order parameter for different values of a. 



a 




2 


h- 


a)/b 


3 


{-a + Va 2 


+ 4&) /2b 


4 


i/\/a 


+ b 



III. RESULTS AND DISCUSSION 

As has been found in previous section, a continuous 
zig-zag transition occurs for parabolic confinement. For 
the case (a ^ 2) it is not possible to define a transi- 
tion between the one-chain and two-chain configuration, 
because the confinement potential docs not contribute to 
the second order term of the total potential energy. Addi- 
tionally we find that the minimum value of the transverse 
frequency is purely imaginary, see Eqs. (|15a[) and (|16a[) . 
and this condition implies that in this case the transition 
for the one-chain to the two-chains configuration is not 
allowed which agrees with previous^ results. We also 
performed Monte-Carlo simulations and found that for 
a > 2 the one-chain configuration is never formed for 
any value of the density and the confinement frequency 
However from similar simulations one can show that for 
a < 2 the one-chain configuration is stable until a criti- 
cal point, beyond which the configuration is changed to 
a single chain containing vacancies due to jumps of indi- 
vidual particles away from the chain axis. 



A. Transition point for a < 2 

For the case of a power-law inter-particle potential 
(k = 0) with parabolic confinement and using dimension- 
less units, it is possible to find an analytical relationship 
between the confinement frequency and the linear den- 
sity as 7/ = u~ 2 A n + 2 ). For this case we show in Fig. @] 
the behavior of the order parameter as a function of the 
linear density for different values of n. Dashed curves 
represent the solution from the Landau theory, Eq. ([5]), 
the full curves are the solution of Eq. (|2"Tj) and they are 
compared with the results of a Monte-Carlo simulation 
for n = 1 and n = 3 (open circles in Fig. From 
these results we notice that there is perfect agreement 
between our calculation and the exact results obtained 
from Monte-Carlo simulations. In the same context Fig. 
[5] shows the variation of the critical density as a function 
of the exponent n as obtained from Landau theory (solid 
and dashed curves for k = and K = 1, respectively) and 
the result from the present work (full and open circles for 
k = and k = 1) the results are shown for different con- 
finement frequencies. Notice that for v = 1 this function 
has a local minimum where the dipole potential exhibits 
the lowest critical density. 

On the other hand, it is also possible to find the value 
of d numerically by fixing one particle at a distance y 
from the one-chain axis in the confinement direction and 



0.0 





















Mr 


/ /" / 


/n=l 


S3 





0.7 0.8 0.9 1.0 1.1 



FIG. 4. (Color online) Displacement from the x-axis as a 
function of the linear density for v = 1. The dashed lines are 
obtained from the Landau theory while the solid lines are the 
solutions of Eq. <|21|) . The open circles represent the results 
of our Monte-Carlo simulations for n = 1 and n = 3 




2 4 6 8 10 

n 



FIG. 5. (Color online) Critical linear density as a function of 
the n exponent of the inter-particle interaction for different 
values of the parabolic confinement frequency. The (solid and 
dashed) lines are the prediction of the Landau theory and the 
(solid and open) circles are found with the present method 
(for k = and k = 1 respectively). 



minimise the energy with respect to the position of the 
other particles. The resulting minimum potential energy 
of the system is shown in Fig. |6] for a Yukawa inter- 
particle potential with k = 1 and v = 1. For ?/ < rj c the 
minimum is found at y = and for i] > r\ c it continu- 
ously shifts to y 7^ 0, which is typical for a second order 
transition. 

From Eq. (|19p we draw the contour plot of v c as a func- 
tion of n and k for several values of the density (when 
making the contour plot we replaced n by a real num- 
ber), which are shown in Fig. [7] We observe a strong 
dependence of the highest value of v c on 77, and therefore 
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0.0 0.2 0.4 0.6 



y 

FIG. 6. Potential energy per particle as a function of the 
position of one of the particles in the confinement direction, 
before (77 = 0.78rj c ), after (77 = 1.34r; c ) and at the transition 
(77 = 77c) for parabolic confinement. We took the parameters 
v = 1, k = 1 and ?i = 1 



the region of frequencies over which the one-chain con- 
figuration exists. For low densities (77 < 1) the one-chain 
organization is dominant for small values of the exponent 
n and gradually this region is extended to higher values 
of n with increasing 77. This result shows that for low 
densities the one-dimensional behavior of the system is a 
better representation for the Coulomb and dipole inter- 
particle potential. For 77 > 1 the one-dimensional region 
of frequencies increasing with increasing n. In all cases 
the critical frequency decreases with increasing re. 

Our previous mean-field theory was derived for modes 
of the linear chain close to the instability point. There- 
fore, it is possible to find the critical point for the afore- 
mentioned instability from Eq. (|27|). In Fig. [8] we plot 
the transition point at which the one-chain structure be- 
comes unstable for a < 2. Above each curve only ip = 
is a solution of Eq. (|27j) . Only for a = 2 the curves 
corresponds to a second order zig-zag transition. Notice 
that the stability region for the single chain configuration 
increases with decreasing a. 

In Fig. [9][a) we show the dispersion relation for the 
normal modes in the case of parabolic confinement for 
different values of k in the three cases, linear regime (dot- 
ted lines) where the system is stable for any value of the 
wavevector close to ko , the zig-zag regime (dashed lines) , 
and in the transition point (solid lines) where the disper- 
sion is linear close to ko. 



B. Case a > 2 

In this case the most simple configuration of the par- 
ticles is restricted to a 2-chains structure, however from 
Monte-Carlo simulations we know that there is a transi- 
tion to a 4-chains structure after some value of the linear 
density. This is shown in Fig. [TUl where we plot the dis- 



tance from the y — axis of the particles as a function of 
the density considering a dipole inter-particle interaction 
for different values of a. In those figures the 2-chains to 
4-chains transition point is marked with a vertical dashed 
line. In our theoretical model we have found from Eqs. 
(|15a[) and (|16a|) that oj\ < as shown in Fig. |9^b) and 
thus the transverse frequency is imaginary and therefore 
the one-chain structure is unstable for any value of the 
density and the confinement strength. 

This is illustrated in more detail in Fig. Qj] where we 
plot the distance of the particles from the x-axis for dif- 
ferent values of 77. Note that our mean-field results from 
Eq. (f27|) agree with the simulation for small values of 
77. Note that for small values the confinement potential 
energy is significantly larger than the inter-particle po- 
tential energy and therefore the fluctuations of the order 
parameter arc smaller. With increasing 77 the interaction 
between the particles start to dominant and all curves 
converge to each other (without crossing) for 77 > 1. 



IV. CONCLUSIONS 

In this work, we studied the critical behavior of a sys- 
tem of particles confined in a 2D channel through a y a 
potential with different functional forms for the inter- 
particle interaction potential. We derived a Ginzburg- 
Landau equation for the system and determine the be- 
havior of the system close to the transition point where 
the single chain configuration becomes unstable. We de- 
termined the order parameter and its dependence on the 
external confinement and the particle density. 

For a = 2 the critical frequency for the zig-zag tran- 
sition is larger than for smaller values of a, which shows 
that the stability of the linear chain configuration is lower 
for parabolic confinement. However for low densities 
(77 < 1) the one-chain configuration is the most stable 
state for a < 2. 

For a > 2 the single chain configuration is unstable for 
any value of the particle density and the strength of the 
confinement potential. We found the distance between 
the two chains as function of the particle density. With 
increasing density a first-order phase transition is found 
to the 4-chains configuration. 

For a < 2 we found analytically no continuous zig-zag 
configuration irrespective of the inter-particle potential. 
The instability of the single chain configuration occurs 
through the expulsion of single particles from the chain 
to y ^ positions. 

The instability point for a = 2 is given by v c = 
^2/(n+2) w j 1 j c j 1 b eC omes an almost linear relation, i.e. 
t) c ~ 77 c for a = 1 and n = 1. 

In a future work we plan to generalize the present anal- 
ysis to the quantum regime, for the special case of elec- 
trons confined by a parabolic potential, i.e. a = 2, n = 1 
and A = 00, such an analysis was presented by J. S. 
Meyer et. ah 30 ' 31 . Subsequently the strongly correlated 
regime which results in Wigner cr ysta l physics in quan- 
tum wires, was addressed in Ref. [32(. Such a quantum 
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FIG. 8. Critical transition point for the one-chain configuration. Above the curves the one chain configuration is the ground 
state. Results are shown for (a) Coulomb potential, (b) dipole potential and (c) Yukawa potential, for different values of the 
exponent of the confinement potential a. 



11 




FIG. 9. (Color online) (a) Dispersion relation lj± as a function of Sk for parabolic confinement (i.e. a — 2) where dashed, solid 
and dotted lines represent the results for v = 0.95u c , v — v c and v = 1.05u c respectively, (b) Square transverse frequency u\ 
as a function of 5k for a > 2. We have considered n = 1, 7? = 1 and different values of k as shown in (a). 




FIG. 10. Distance of the particles from the y-axis as a function of the density, for different values of a and v = 1. The 
inter-particle interaction is a dipole potential. The symbols represent the results from Monte-Carlo simulations and the dashed 
vertical lines indicate the transition point from the 2-chains to the 4-chains configuration. 



analysis will address the effect of quantum statistics of 
the particles and the effect of quantum fluctuations on 
the zig-zag transition. 
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